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Abstract 

' We develop an experimental algorithm for the exact solving of the 

maximum independent set problem. The algorithm consecutively finds 
the maximal independent sets of vertices in an arbitrary undirected graph 
such that the next such set contains more elements than the preceding one. 
For this purpose, we use a technique, developed by Ford and Fulkerson for 
^ ' the finite partially ordered sets, in particular, their method for partition 

| of a poset into the minimum number of chains with finding the maximum 

antichain. In the process of solving, a special digraph is constructed, and 
a conjecture is formulated concerning properties of such digraph. This 
allows to offer of the solution algorithm. Its theoretical estimation of 
running time equals to is 0(n 8 ), where n is the number of graph vertices. 
' The offered algorithm was tested by a program on random graphs. The 

, testing the confirms correctness of the algorithm. 
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1 Introduction 



A set of all undirected n-vertex graphs without loops and multiple edges is 
denoted by L n . 

Let there be a graph G = (V, T) E L n , where V is the set of graph vertices, 
and r is a mapping from V to V. Any subgraph Qi = (Vi,Ti) of G is called a 
clique induced by V\, if 

T 1 (v) = V 1 \{v} 

for all v € Vi. In a special case, when Card{V\) = 1, the one- vertex subgraph 
Qi = (Vi,Ti) is called a single-vertex clique. 

A clique Q\ is called maximal if any vertex v G V cannot be attached to it 
so that the new vertex set also has formed a clique of the graph G. A clique Q 
is called maximum if the graph G has not a clique of the greater size than Q. 
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It be required to find the maximum clique of a graph G G L n . The known 
problem of finding the maximum clique has formulated [3]. 
A vertex set U C V of G is called independent if 

UC\T(U) = 0. 

An independent set U of graph vertices is called maximal (MIS) if 

U U T(U) = V. 

A MIS U is called maximum (MMIS) if Card(U) > Card(U) for any MIS 
U of G. 

It be required to find the MMIS of a graph G. Again, we have formulated 
well- know the maximum independent set problem (MISP). 

Both of the formulated above problems are NP-complete [31 Q] . They are 
closely connected with each other: the solution of one of them attracts the 
solution another. 

A graph G is called complement to the graph G if it has the same vertex set, 
and edges join two vertices of the graph G iff these vertices are non-adjacent in 
G. 

It is not difficult to see that any clique of G corresponds to the independent 
set of graph vertices in G, and conversely. Therefore, the finding the maximum 
independent set in one graph is equivalent to the finding the maximum clique in 
its complement graph. In the given work, we examine the maximum independent 
set problem. 

The intention of the given paper is to design a polynomial-time algorithm 
for the exact solving of the maximum independent set problem in an arbitrary 
undirected graph. For this purpose, we use a technique, developed by Ford 
and Fulkerson (see [2] and Appendix A) for the finite partially ordered sets, in 
particular, their method of partition poset into the minimum number of chains 
with finding the maximum antichain. In the process of the solving, a special 
digraph is constructed, and a conjecture is formulated concerning properties 
such digraph. This allows to offer a solution algorithm, a theoretical estimation 
of running time for which equals to is 0(n 8 ), where n is the number of graph 
vertices. The offered algorithm was tested by a program. The testing confirms 
the correctness of the algorithm. 

Notice that the author did not aspire to creation of the most optimum and 
fast algorithm. In work [8j, Tarjan has presented a table, which shows a dynam- 
ics of perfecting a complexity evaluation for solving some problems. Whence 
one can make a conclusion that after the appearance of an initial algorithm for 
solving certain problem, its improvement is found quickly. It is presented that 
proposed algorithm also can be improved hereinafter. 
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2 The Basic Definitions 



Let there is a graph G = (V, T) 6 L„. We will partition the vertex set V into 
subsets 

V°, V 1 , . . . , V m (1) 

in such a way that a subset V k (k = 0,1, ... ,m) is a MIS of a subgraph Gk 
= (V \ (V° U V 1 U • ■ ■ U V^ 1 ), r fe ) = (V k U • ■ ■ U V m , r fe ). Clearly, G = 

(^u--ur, r ) = g. 

By the given undirected graph G and the partition ([T]) we can construct a 
digraph G(V°) = (V, T) in the following way. If an edge of G joins a vertex 
Vi G y' 01 with a vertex Vj G V k2 then this edge is replaced by an arc (vi,Vj) 
when k\ < k2- The vertex Vi is called the tail of (vi,Vj), and the vertex Vj is 
called the head of this arc. 

As the result we have an acyclic digraph G(V°) = (V,T). The set V° is 
called initiating. 

In general case we can construct a set D(G) of different acyclic digraphs as it 
was indicated above. Each digraph of D(G) corresponds to the graph G G L n . 
Further we will consider only digraphs of D(G). 

The maximum length p(v) of a directed path, connecting a vertex v G V 
with some vertex of the initiating set V°, is called the rank of v. The set of all 
graph vertices having the same rank p(v) = k is called the k-th layer of G(V°) 
and designated as V k . 

To apply the partially ordered set technique, each digraph G(V°) is assigned 
to a transitive closure graph (TCG) <3 t (V°) = (V,f t ) [§]. As the digraph 
G(V°) is acyclic and loopless, its transitive closure G t (V°) is a graph of a strict 
partial order (V, >). Further, we will not distinguish the transitive closure graph 
Gt(V°) and partially ordered set (poset) (V, >). Therefore, we will consider, for 
example, antichains of the TCG Gt(V°). 

There exists an efficient algorithm to construct the TCG. Its running time 
is equal to 0(n 3 ) (see, for example, [7J [B]). 

An arc (vi,Vj) oiGt(V°) will be called essential if there exists the arc (vi,Vj) 
of the digraph G(V°). Otherwise, the arc (vi,vj) will be called fictitious. An 
essential arc is also designated as Vi > Vj , and a fictitious arc is also designated as 
Vi 3> Vj . Obviously that any fictitious arc of Gt (V°) determines two independent 
vertices of the digraph G(V°). 

Let there is a poset (A, >). 

If a > b or b > a, the elements a and b of A are called comparable. If a ^ b 
and b a, such pair of elements is called incomparable. 

If A\ C A and each pair of elements of A\ is comparable, we shall say that 
A\ determines a chain of (A, >). If A\ G A and each pair of elements of A\ is 
incomparable, we shall say that A\ be an antichain of (A, >). The antichain Ai 
is the maximum in (A, >), if Card(Ai) > Card{A*) for any antichain A* G A 
in(A,>). 
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We say that poset (A, >) is partitioned into chains Ai, . . . , A p , if each Ai 
(Ai ^ 0, i = be a chain, 

U A: = ^, 

and ^ n Aj = 0, when i 7^ j (i,j 6 {1, . . . ,p})- 

The partition of the poset (A, >) into chains is called minimum, if it has the 
minimum number of elements p in comparison with other partitions of (A, >) 
into chains. Such partition also is called minimum chain partition (MCP) of 
poset (A, >). 

As G t (V°) is a graph of strict partial order, we can find MCP V = {Si, . . . , S p }. 
In common case, this partition is ambiguous. 



v 5 v 6 v 5 v 6 




Vl V 2 Vl v 2 

(a) (b) 
Figure 1: Different MCPs of the transitive closure graph 

Different MCPs of the transitive closure graph is shown on the Fig. [T] (a) 
and (b). The digraph arcs, belonging chains of MCP, are represented by thick 
lines. Here and hereinafter, we suppose that orientation of arcs of the digraph 
is from below to upwards. 

Let V(Sq) (q = l,p) be the set of vertices, belonging to the chain S q . If 
vertices Vi, vj 6 V(S q ) are endpoints for a fictitious arc w,; 3> Vj then the vertex 
Vj is called marked. The set of all marked vertices of Gt(V°) is determined 
by the found MCP V and it differs for different MCPs. The set of all marked 
vertices of G t (V°) is designated as B(V). 

For example, for the MCP Vi represented in Fig. [T] (a), we have B(Vi) = 
{V5, v 6 }, and for the MCP V2 represented in Fig. Q] (b), we have B(V2) = 0- 

Lemma 1 Let G(V°) be a digraph, and Gt(V ) be its transitive closure graph. 
If B(V) = 0, where V is an MCP of Gt(V°), then the maximum antichainU is 
the MMIS of the graph G G L n , and the MCP determines the minimum clique 
partition of G. 

If conditions of Lemma [1] are satisfied then each chain S q G V is a clique of 
G(V°). Therefore, the MCP V is the minimum clique partition. 
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On the other hand, vertices of the maximum antichain U of Gt(V°) belong 
to distinct cliques, that is, the number of vertices in the MMIS is equal to the 
number of vertices in the set U. Q.E.D. 

Notice that Lemma[T]may be satisfied if the digraph G(V°) has no transitive 
orientation. 

It is obvious that any antichain U C V of Gt(V°) determines an independent 
vertex set of the digraph G(V°), and an independent vertex set U C V of G(V°) 
determines an antichain of Gt(V°) if and only if no two vertices of U belong to 
the same directed chain of G(V°). 

3 A vertex-saturated digraph 

Let there is an acyclic digraph G(V°) — (V, T). 

Further, let W C V be some independent vertex set. For the digraph G(V°) 
we define an unary operation cutting o~w(G(V )). This operation consists of 
reorientation of all arcs of G(V°) incoming into vertices of the set W. It is easy 
to see that the result of this operation is also a digraph G(Y°), where 

r° = (v°\f- 1 (w))uw. 

Here, T _1 is a maping, inverse to T. 

Theorem 1 Let there be a digraph G(V°) S D(G) and W be some independent 
vertex set. Then a digraph G(Y ) = o~w(G(V )) is also acyclic and G(Y°) S 
D{G). 

Indeed, since the digraph G(V°) = (V,T) is acyclic then any its part is also 
an acyclic digraph. Thus, a directed subgraph Gi = (V \ W, Ti) is acyclic. 

Obviously, G\ C G(Y°). Attach the independent vertex set W to the sub- 
graph G\. Join each vertex x 6 W with a vertex y of G± by the arc (x, y) if and 
only if there exists the arc (y,x) of the digraph G(V°). It is evident that the 
resulting digraph G(Y°) is also acyclic. 

At last, we have G(Y°) 6 D(G) since any reorientation of arcs of the digraph 
G(V°) does not change independence relation of its vertices. Q.E.D. 

Let there be a digraph G(V°) and its transitive closure graph G t (V°). 

We can find the MCP V of the graph G t (V°) constructing the maximum 
antichain IA simultaneously as it is described in Appendix \^ 

In general case we can find some distinct maximum antichains. 

We will say that an antichain U\ precedes an antichain Ui in the graph 
Gt(V°) and designate it as U\ ~< U2 if for all vertices x € Ui \ U2 there is a 
vertex y 6 U2 \ U\ such that x < y. 
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By means of Ford and Fulkerson's methodology such a maximum antichain 
U of the TCG Gt(V°) can be found that precedes any other maximum antichain 
Ui, i.e. U\ -n U for any antichain U\ of the TCG Gt(V°). This antichain of the 
graph Gt(V°) we will call general. 

In addition to the general antichain, we may find other maximum antichains 
of G t (V°) if they exist. So for any vertex v G V of the TCG G t (V°), it is possible 
to find a maximum antichain U{v) such that v G U(v). Technically, to find the 
antichain U(v), it is sufficient, in the adjacent matrix of Gt(V°) containing the 
maximum number of units in allowable cells (and the marks are appointed by 
the Ford-Fulkerson's algorithm), to add the mark (*) to the existing marks for 
a row, corresponding to vertex v, and to execute a cycle of appointment of 
marks. In this case, in the first step of appointment of marks, all columns are 
marked, which contain the admissible cells (incluging a chosen cell). Clearly, 
the antchain U(v) will be general for the vertex v, that is, any other antichain, 
containing vertex v, will precede this antichain. 



a.2 aj, 



«5 




Figure 2: The partially ordered set 

Notice that in general case it is not true that for any vertex v G V of the 
TCG Gt(V°) there exists a maximum antichain U(v). For example, consider 
the digraph of the partially ordered set, shown in Fig. [2J It is easy to see that 
there are no maximum antichains U(a\) and U{ai) in it. 

A directed subgraph G(V k ) (k — 0, m — 1) of the digraph G(V°) we will 
call saturated with respect to the initiating set V k , if, in its transitive closure 
graph Gt(V k ), any maximum antichain U(v) C V, when it exists, is a MIS of the 
subgraph G{V k ) and satisfies the relation: Card(U{v)) — Card(V k ). Evidently, 
the subgraph G(V m ) has no arcs, and therefore it is saturated with respect to 
its initiating set always. 

A digraph G(V°) is called vertex-saturated (VS-digraph) if any of its directed 
subgraphs G(V k ) (k — 0, m), induced by the layer V k , is saturated with respect 
to the initiating set V k . 

Notice that digraph, represented in Fig. [lj is not vertex-saturated. 

Let there be some digraph G(V ). To construct a VS-digraph, we use the 
following algorithm. 

The algorithm VS. 



G 



Step 1. Put k := and a := false. 

Step 2. Find the transitive closure graph G t (V k ). 

Step 3. Construct an MCP of G t {V k ). 

Step 4- Find the maximum antichain of G t (V k ) for each vertex v G V k U • • • U 
V m . 

Step 5. Check whether each of the found maximum antichains U(v) of the graph 
G t {V k ) is a MIS of the digraph G(V k ) and Card(U{v)) = Card(V k ). If 
it is true, finish the design of the digraph G(V k ) saturated with respect 
to the initiating set V k . Go to Step 6. 

Otherwise complete the found antichain U(v) (when it is necessary) to 
a MIS, put W :— U(v) and construct a new acyclic digraph G(V k ) by 
the cutting operation aw(G(V )), put a :— true and construct the new 
digraph G(V°). Return to Step 2. 

Step 6. Compute k := k + 1. If k < m then distinguish the transitive closure 
graph G t {V k ) from G t (V fc_1 ), keeping all chains of the MCP of G t (V k - 1 ) 
which are incident to vertices of the new graph. Go to Step 4. 

If k = m and a = true, go to Step 1. If k = m and a — false, go to Step 
7. 

Step 7. Finish of this algorithm. A VS-digraph is constructed. 

We will show that the algorithm VS constructs a vertex-saturated digraph. 

Theorem 2 Let G(V ) be a digraph constructed by the algorithm VS. Next, let 
G(V k ) = (Yfc,r) (fc = 0, m) be a directed subgraph induced by the layer V k . 
Then each antichain U C Yfe \ V k of the graph G t {V k ) — {Yk,Tt) obeys the 
following relation 

Card{V k n T^U) > Card(U) (2) 

Assume that, in the directed subgraph G(V k ), an antichain U C Yk \ V k of 
G t {V k ) = (Yfc,r t ) will be found such that 

Card(V k n Tt l $l)) < Card(U). 

Construct a set U* = (V k \ r t _1 (W)) UW. Clearly, the set U* is independent 
in the directed subdigraph G(V k ) and CardiU*) > Card(V k ). 

Since the set V is an antichain of the TCG G t {V k ) = (Y k ,T t ) then the set 
V k \ r t _1 (W) C V k is an antichain of this graph. The set U C Y k \ V k is an 
antichain of Gt(V k ) by conditions of Theorem [2] 

Obviously, U H T t (V k \ r t _1 (W)) = 0, therefore the set W is an antichain of 
G t (V k ). 

We have obtained a contradiction since the directed subgraph G(V k ) is con- 
structed by the algorithm VS and the antichain Li* could be discovered in Step 
5. This proves the validity of Theorem H Q.E.D. 
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Corollary 1 The digraph G(V°), constructed by the algorithm VS, is vertex- 
saturated. 

Theorem 3 VS-digraph can be constructed in time 0(n 5 ). 

One completion of the steps 1-6 requires 0(n|) time units, where nk is 
the size of the vertex set of G(V k ). Assuming that for each completion of these 
steps the maximum antichain increases at one vertex, we obtain the design time 
of a digraph G(V k ) vertex-saturated with respect to the initiating set, is equal 

to oK)- 

Hence, the design time of a correctly constructed vertex-saturated digraph 
G(V°) is equal to: 

Y,0(ni)=0(n 5 ). 

Q.E.D. 

Theorem 4 Let a digraph G(V ) be vertex- saturated. Then there exists an 
MCP V of the graph Gt(V°) such that its chains contain only essential arcs. 

Let a digraph G(V°) is constructed by the algorithm VS. By Theorem^ 
each bipartite digraph G(V k , V k+1 ) (k — 0, m) of this digraph satisfies the Hall's 
theorem and, hence, has a matching that saturates each vertex of the set V fc+1 ). 
Q.E.D. 

Corollary 2 Let G{V°) be a VS-digraph. Then each chain of an MPP V of 
Gt(V°) is begun by some vertex v ofV . 

Due to this result we may use the adjacent matrix of G(V°) as a working 
table for determination MCP of the TCG G t (V°) of a VS-digraph. Thus, we 
will suppose that chains of MPP of the TCG of a VS-digraph are found using 
the adjacent matrix of this digraph. That is, we choose only essential arcs 
to construct each new MPP! 

Certainly, we use the adjacent matrix of the transitive closure graph to search 
for the maximum antichains U (v) of such TCG. 

The instance of constructing vertex-saturated digraph is shown in Appendix 

m 

4 An algorithm for finding MMIS of a graph 

Let a saturated digraph G(V°) is constructed, which has a MMIS U such that 
Card(U) > Card(V°). In this case, at least one of the chains of TCG of the VS- 
digraph G(V°) contains a fictitious arc, whose endpoints belong to the MMIS. 

Let, further, a some fictitious arc Vi ^> Vj is found in the TCG Gt(V°). 
We shall remove the vertices Vi,Vj from the digraph Gt(V°) and all vertices, 
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which are adjacent with them. As a result, we shall obtain a digraph Gi(V 1 ) = 
(Vi,Ti), where 

V 1 =V\({v i ,v j }UT(v i )UT(v j )), 

v? - v° \ (f-\ Vi ) uf- 1 ^)). fa =f nv x . 

HereT(w) = f(v) U f- x (u). 

For the digraph Gi(V]°), we shall use the procedure of constructing a VS- 
digraph by the algorithm VS. As a result, we shall obtain a digraph G(Z°), 
which shall call induced by removing the fictitious arc Vi 3> Vj . 

An algorithm for finding a MMIS of a digraph G(V°) is constructed on the 
supposition that the following conjecture is true. 

Conjecture 1 Let a saturated digraph G(V°) has an independent set U C V 
such that Card(U) > Card(V a ). Then it will be found a fictitious arc Vi 3> 
Vj such that in the digraph G(Z°), induced by removing this arc, the relation 
Card(Z°) > Card(V a ) - 1 is satisfied. 

The worded conjecture allows to formulate a solution algorithm for finding 
a MMIS of a graph G € L n . Input of the algorithm is an undirected graph 
G e L n . Output of the algorithm is the MMIS. 

An algorithm for finding a MMIS. 

Step 1. Execute an initial orientation of the graph edges so to get an acyclic 
digraph G(V°). 

Step 2. Execute the algorithm VS for the digraph G(V°). 

Step 3. In TCG of the VS-digraph to find an unmarked fictitious arc Vi 3> Vj. 
Mark the found fictitious arc as considered. If all fictitious arcs are marked, 
go to the Step 7. 

Step 4- Remove vertices Uj, Vj as well as all adjacent with them vertices from the 

saturated digraph G(V°). As a result, a digraph G 1 (V r 1 °) will be obtained. 
Step 5. Execute the algorithm VS for the digraph G\ (V]°). As a result, a digraph 

G(Z°) will be obtained. 
Step 6. If Card(Z Q ) > Card{V°) - 1, construct a set W = Z° U {v^Vj} and 

execute the cutting operation aw(G(V a )) in the saturated digraph G(V°). 

Go to Step 2. Otherwise go back to Step 3. 

Step 7. Put a MMIS U — V°. 

Theorem 5 If the conjecture [7] is true then the stated algorithm finds a MMIS 
of the graph G S L n . 

It is obviously. Q.E.D. 

Theorem 6 The running time of the algorithm of finding a MMIS eguals to 
0(n s ). 
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Indeed, single executing of the Steps 3-6 requires 0(n\) of time units, 
where n\ is the number of vertices in the digraph G(Z ), induced by removing 
a fictitious arc. Since total number of fictitious arcs is 0(n 2 ), in worse case, 
for executing the Steps 3 - 6 is required 0(n 7 ) of time units. If suppose that 
after executing these steps, the found independent set V° will be increased to 
the unit, the total running time of steps 2-6 equals to 0(n 8 ). Q.E.D. 

5 Conclusion 

The pascal-programs were written for the proposed algorithm. Long testing 
the program for random graphs has shown that the algorithm runs stably and 
correctly. 

Of course, the offered algorithm is not competitive in practice because of high 
degree of polynomial estimation of the running time. However, the algorithm has 
important theoretical significance since there is a good probability to prove that 
for NP-complete problems possible to construct a polynomial-time algorithm. 

Author thanks Guenter Stertenbrin!0 for the help in testing the programs. 
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Appendix 



A Partially Ordered Sets 



We recall some conceptions of Set Theory. 

Relations between two objects are called binary. A binary relation R may 
be represented by a listing of object pairs, which are in the relation R: 

(oi,&i), • ■ • , (a m ,b m ). 

If (a, b) G R, then this fact we also denote by aRb. When (a,b) $ R, then we 
will write aRb. If a S A and b ^ A for all aRb, then i? is called a relation on 
the set A. Further, we will consider only relations on the finite set A. 

If A is a finite set and R is a relation on A, we can represent i? as a digraph 
G. Each element of A is assigned to a vertex of G, and the vertex is joined 
with a vertex a, by the arc (a.i, a,-) if and only if aiRaj. 

A relation R is reflexive if ai?a for every a 6 A A relation R is irreflexive 
if aRa for every a 6 A. A relation i? is symmetric if whenever ai?6, then bRa. 
A relation i? is antisymmetric if whenever ai?6 and bRa, then a = b. A relation 
i? is transitive if whenever ai?6 and bRc, then ai?c. 

A binary relation R is called a partial order if i? is antisymmetric, and 
transitive. The set A together with the partial order R are called a partially 
ordered set. We will denote this partially ordered set by (A, R) or (A,>). If 
a relation R is irreflexive, then such partial order is called strict. A strictly 
ordered set is written by (A, >). 
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Figure 3: The adjacent matrixes 

In Fig. [21 the acyclic graph represents the partial order, induced by the 
binary relation 

R = {(ax, <z 2 ), (01,03), (01, 05), (a 4 , a 2 ), (a 5 , a 2 )}. 

Here and throughout, we assume that the orientation of arcs of a digraph on a 
drawing is from below upwards. 

Dilworth's famous theorem establishes the relationship between a MCP and 
the maximum antichain of (A, >) [HH]. 
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Theorem 7 (Dilworth R.P.) Let (A, >) be a finite partially ordered set. The 
minimum number of disjoint chains, which the set (A, >) can be partitioned on, 
equals to the capacity of the maximum antichain in (A, >). 

There is an efficient algorithm for the partitioning a finite partially ordered 
set into the minimum number of chains and for finding the maximum antichain, 
elaborated by L. R. Ford and D. R. Fulkerson ,2 ]. In essence, this algorithm 
finds the maximum matching in a bipartite graph G* = (X, Y, T*). If a partially 
ordered set has n elements, then this graph contains 2n vertices and Card(X) = 
Card(Y) — n. An edge (xi,yj) joins two vertices x\ G X and yj S Y if and 
only if the corresponding elements a%, a-i £ A of (A, >) are comparable. 

In manual computations, we will use an adjacent matrix M of G* as a 
working table. Units of M determine its admissible cells. Two cells of M are 
called independent if they are located in distinct rows and distinct columns of 
M. To find the maximum matching, we will have to find the maximum number 
of admissible independent cells of M. 

The algorithm for partitioning a partially ordered set into the minimum 
number of chains consists of two stages: 

• Construct an initial partition of the partially ordered set into chains; 

• Improve the existing partition if it is possible. 

To avoid ambiguity, we always look through rows and columns of M uni- 
formly: from top to bottom in columns and from left to right in rows. 
To obtain an initial partition, we may use the following procedure. 

Step 0. Put N = n, where n is the number of rows M, i = 1. 

Step 1. If N = 0, then complete the calculations as an initial partition is found. 

Step 2. In i-th row of M find the first on the order admissible cell, whose ap- 
propriate column is not marked. If such cell is not found, put i := i + 1, 
N := N — 1 and go to Step 1. Otherwise, remember the found cell (i,j), 
mark a column j, calculate i := i + 1 and go to Step 1. 

Fig. [3] (a) shows the adjacent matrix for the strictly ordered set, represented 
in Fig. [5) The chosen cells of the initial partition are indicated by a circle. 

To find a MCP and the maximum antichain of a set (A, >), we will make use 
of the Ford-Fulkerson's algorithm .2] . The algorithm begins to work after termi- 
nation of the previous procedure, that is, when there exists an initial partition 
of the ordered set into chains. 

Step 1. Mark rows of M that do not contain the chosen cells, by the symbol 
(*)• 

Step 2. Look through the newly marked rows of M and find all unchosen cells 
in each row. Mark all unmarked columns of M that correspond with such 
cells by an index of the row. 
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Step 3. Look through the newly marked columns. If an examined column con- 
tains a chosen cell (that is, the cell is enclosed within a circle), then mark 
the row containing the chosen cell by an index of the examined column. If 
the column does not contain a chosen cell, go to Step 4. If it is impossible 
to mark new rows, then go to Step 5. 

Step 4- The essence of the given step is the procedure of constructing a new 
collection of independent cells, each having one more cell than the former 
collection. At each stage of this procedure, except for the final step, we 
pick a new admissible cell of M and delete the "old" one. Increasing 
the total number of chosen cells happens as follows. In the found j-th 
column, choose a new cell in a row m(j), where m(j) is a mark of the 
current column. Let we already have chosen the cell (i,j), which marks 
m(i) and m(J) correspond to, where m(i) is a mark the i-th row. If 
m(j) = (*), the procedure of constructing a new collection of independent 
cells is completed. Delete all marks of rows and columns, and go to Step 
2. Otherwise, delete the cell (i,m(i)) and choose a cell (m(m(i)), m(i)). 
Put i = m(m(i)), j = m(i) and repeat the process described above. 

Step 5. Find the maximum antichain U — U r \U c , where U r is a set of marked 
rows, and U c is a set of marked columns. Terminate the calculations. The 
found admissible cells determine arcs forming chains of the MCP. 

Fig. [3] (b) shows the picked cells of the optimal partition for the partially 
ordered set, represented in Fig. [21 In this case, the MCP consists of the 
chains A\ — {ai,a 3 }, A 2 = {04,02}, and A3 = {a 5 }. We also have a set 
U r = {02,03,04,05} of marked rows, and a set U c — {02} of marked columns. 
Consequently, the maximum antichain 11 is equal to 

U = {02,03,04, a 5 } \ {a 2 } = {a 3 , a 4 , a 5 }. 

Ford-Fulkerson's methodology, described above, for finding antichains may 
be easily adapted to any algorithm of finding the maximum matching in a bi- 
partite graph, for example, Hopcroft-Karp's algorithm [J], [7j, or flow algorithm 
[5 . Therefore, we assume that the running-time of a MCP construction is equal 
to 0(n 5 / 2 ). 

B An instance of constructing a vertex-saturated 
digraph 

Consider an instance of constructing a vertex-saturated digraph (VS-digraph). 

Fig. [5] shows a digraph G(X°), obtained from an initial undirected graph G 
as it was described in part O Recall that the orientation of arcs of the digraph 
in figures is from below upwards. 

Construct a VS-digraph. 
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The adjacent matrix of the TCG Gt(V°) is shown in Fig. [5j The fictitious 
arcs of this graph are represented by the letter f. Arcs, belonging to the MCP 
of Gt(V°), are put into circles. These arcs are shown in Fig. 0]by thick lines. 

First of all, notice that the initiating set V° = {v\, v 2 , V3, V4} is a MIS of 
G(V°). 

Find the general antichain hi (see Fig. [5] (a)). The set of marked rows is 
U r = {5,7,8,10}, and the set of marked columns is empty, that is, U c = 0. 
Therefore: 

U = U r \U c = {5, 7, 8, 10} \ = {5, 7, 8, 10}. 

This antichain is a MIS of the digraph, and Card(U) = Card(V°). 
Find the maximum antichains U(v) for graph vertices. 
Clearly, U(v 5 ) = U{v 7 ) = U(v s ) = U{v w ) = U. 
To find U(vi), mark the first row of the adjacent matrix of Gt(X°). 
Marking the first row in Fig. [5](b), we have U r = {1, 3, 4, 5, 6, 7, 8.9.10}, and 
U c = {5,7,8,9,10}. Consequently, 

U(vi) =U r \U c = {1,3,4,6}. 

This antichain is a MIS of the digraph, and Card(U(vi)) = Card(X°). 
Similarly, marking the second row in Fig. [5] (c), we have 

U(v 2 ) = {1, 2, 3, 4, 5, 6, 7, 8, 9, 10} \ {5, 6, 7, 8, 9, 10} = {1, 2, 3, 4}. 

This maximum antichain is also a MIS of the digraph, and Card(U(v2)) = 
Card(V°). 

Similarly, we obtain the maximum antichains 

U(v 3 ) = {3,4,5,8}, 
U{v A ) = {4,5,8,9}, 

and 

U(v 6 ) = {3,4,5,6} 




vi v 2 v 3 V4 



Figure 4: A digraph 
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Figure 5: Finding antichains 
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Figure 6: The directed subgraph induced by the layer V 1 
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Figure 7: The new directed subgraph 



from Fig. 0(d), (e), and (f) correspondingly. Each of these maximum antichains 
is a MIS of the digraph G(V°), and they have the number of elements equal to 
V°. 

Thus, the digraph G(V°) is saturated with respect to the initiating set V°. 

Now we examine directed subgraphs induced by layers of G(V°). 

Consider directed subgraph G^ 1 ) induced by layer V 1 = {v^,vq,v 7 }. This 
subgraph is represented in Fig. [6] (a). Notice that the adjacent matrix of the 
TCG Gt{V r ) can be obtained from the adjacent matrix of Gt(V°) directly. 
Obviously, the MCP of GtiV 1 ) is a part of the MCP of G t (V°). The adjacent 
matrix of GtiV 1 ) is shown in Fig. [6] (b). The general antichain U\ of GtiV 1 ) 
equals 

Ul = {«5,« 7 ,«8,Uio}. 

This antichain is a MIS of G(V r ); however, Card{Ux) > CardiV 1 ). Con- 
sequently, the directed subgraph GCV 1 ) is not vertex-saturated with respect to 
its initiating set. 

Therefore, we assume W — {v§, v 7 , vs, Vio} and reorientate all arcs of G{X 1 ) 
incoming to the vertices of W. We obtain a new directed subgraph represented in 



17 



Fig. [7] (a). Clearly, this subgraph has a new initiating set V 1 — {v$, V7, v$, vio}. 

Examining as above, we find that the new directed subgraph G^V 1 ) is vertex- 
saturated with respect to its initiating set V 1 . 

Thus, we may construct a new digraph G(V°). The adjacent matrix of this 
digraph can be obtained from the adjacent matrix of the initial digraph if the 
corresponding part of it is replaced by the adjacent matrix of GfV 1 ). 

Similarly, we determine that a directed subgraph G(V 2 ), where V 2 — {6}, 
is vertex-saturated with respect to its initiating set V 2 . 
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Figure 8: The VS-digraph 



At last, we may make sure that the new digraph G(V°) is a VS-digraph since 
each of its directed subgraphs is vertex-saturated with respect to its initiating 
set. This digraph is represented in Fig. [8] (a). The adjacent matrix of the 
transitive closure graph G t (V°) and its MCP together are represented in Fig. [8] 
(b). 
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